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The AGILE satellite, operating since mid 2007 and primarily devoted to high-energy astrophysics, 
is one of the only three currently operating space instruments capable of detecting Terrestrial 
Gamma-Ray Flashes (TGFs), together with RHESSI and Fermi-GBM.. Thanks to the AGILE 
Mini-Calorimeter instrument energy range extended up to lOOMeV and its flexible trigger logic on 
sub-millisecond time scales, AGILE is detecting more than 10 TGFs/month, adding a wealth of 
observations which pose severe constrains on production models. The main AGILE discoveries in 
TGF science during two and a half years of observations are the following: 1) the TGF spectrum 
extends well above 40 MeV, 2) the high energy tail of the TGF spectrum is harder than expected 
and cannot be easily explained by previous theoretical models, 3) TGFs can be localized from space 
using high-energy photons detected by the AGILE gamma-ray imaging detector. In this presenta- 
tion we will describe the characteristics of the 2.5-years AGILE TGF sample, focusing on the recent 
results concerning the TGF high-energy spectral characteristics. 
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INTRODUCTION 



Terrestrial Gamma-Ray Flashes (TGFs) are very 
short (lasting up to a few milliseconds) bursts of 
high-energy photons above 100 keV, first detected by 
the BATSE instrument on board the Compton Ob- 
servatory [![, and later extensively observed by the 
RHESSI satellite @, i]. TGFs have been associated 
with strong thunderstorms mostly concentrated in the 
Earth's equatorial and tropical regions, at a typical al- 
titude of 15-20 km Q. TGFs are widely believed to be 
produced by Brcmsstrahlung in the atmospheric lay- 
ers by a population of runaway electrons accelerated 
to relativistic energies by strong electric fields inside or 
above thunderclouds. The secondaries generated dur- 
ing the acceleration process can be accelerated as well 
driving an avalanche multiplication Q , commonly re- 
ferred to as Relativistic Runaway Electron Avalanche 
(RREA) , possibly further enhanced by means of a rel- 
ativistic feedback mechanism @, 0j- Although being 
widely accepted as the underlying physical process in 
TGF production, the RREA mechanism alone is not 
sufficient to explain the rich phenomenology of TGFs, 
especially the observed fluence, and there is no consen- 
sus yet on the underlying physical conditions, produc- 
tion sites, radiation efficiencies and maximal energies. 



Recently, the AGILE satellite added a wealth of 
spectral and geographical data on TGFs and estab- 
lished itself as a major player in TGF observation, 
together with the RHESSI || and Fermi-GBM de- 
tectors j8|. 

AGILE 9] is a mission of the Italian Space Agency 
(ASI) dedicated to astrophysics in the gamma-ray en- 
ergy range 30 MeV - 30 GeV, with a monitor in the 
X-ray band 18 keV - 60 keV [l(J, operating since 
April 2007 in a low inclination (2.5°) equatorial Low- 
Earth Orbit at 540 km altitude. The AGILE Gamma- 
Ray Imaging Detector (GRID) is a pair-tracking tele- 
scope based on a tungsten-silicon tracker [Tl|, sen- 
sitive in the 30 MeV - 30 GeV energy range. The 
imaging principle is based on the reconstruction of 
the tracks left in the silicon detection planes by the 
electron-positron pairs produced by the primary pho- 
ton converting mainly in the tracker tungsten planes. 
A Mini-Calorimeter (MCAL) [Hj], based on CsI(Tl) 
scintillating bars for the detection of gamma-rays in 
the range 300 keV - 100 MeV, and a plastic anti- 
coincidence detector [l3[ complete the high-energy in- 
strument. MCAL can work also as an independent 
gamma-ray transient detector with a dedicated trig- 
ger logic acting on several time scales spanning four 
orders of magnitude between 290 /is and 8 seconds 
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GIG!. 

Thanks to its flexible trigger logic on sub- 
millisecond time scales [3], MCAL proved to be a 
very efficient instrument for TGF detection. In this 
paper we will review the AGILE achievements in TGF 
science, focusing on our results on TGF high-energy 
spectral characteristics and on the correlation with 
global lightning activity. 

II. AGILE RESULTS IN TGF SCIENCE 

The average MCAL detection rate is 
~ 10 TGFs/month, with the current severe se- 
lection criteria based on hardness ratio and fluence 
[ltij . When a trigger is issued by the onboard logic, 
data are collected on an event-by-event basis so that, 
for every photon, energy, timing with 2/xs accuracy 
and topological informations (i.e. the fired detector 
address) are saved and sent to ground. For a trigger 
to be classified as a valid TGF candidate, at least 
10 photons and a hardness ratio HR > 0.5 are re- 
quired, where HR is defined as the ratio between the 
number of counts with energy greater than 1.4 MeV 
and the number of counts with energy lower than 
1.4 MeV. There is evidence that relaxed selection 
criteria may be applied, resulting in the detection of 
a fainter/softer population which can increase the 
detection rate more than 50% (Marisaldi, M. et at, 
in preparation). 

The main AGILE discoveries in TGF science during 
two and a half years of observations are the following: 

- the TGF spectrum extends at least up to 
40 MeV [Hi (well above the previous 20 MeV 
limit set by RHESSI @); 

- the high energy tail of the TGF spectrum is 
harder than expected and cannot be easily ex- 
plained by previous theoretical models [17j : 

- TGF can be localized from space using high- 
energy photons detected by the AGILE silicon 
tracker [lj] 

- TGFs are not a random sub-sample of global 
lightning activity as detected from space (l9j . 
Moreover, significant regional differences exist, 
both in the degree of correlation and in the 
TGF/flash ratio. 

A. TGFs high-energy spectrum 

The cumulative spectrum of the first 34 high- 
confidence TGFs detected by AGILE shows a good 
agreement with the RHESSI cumulative spectral 
shape below 20 MeV, when the detector response ma- 
trices for both instruments are taken into account, but 



demonstrates also that events of energy up to 40 MeV 
can be observed [l6j]. These early observations sub- 
stantially doubled the previously known TGF energy 
range set by RHESSI as later confirmed by the 
observations of Fermi-GBM 0] . 

The high energy spectral shape of TGFs became 
even more puzzling when increased statistics was 
achieved. Fig. [1] shows the background-subtracted 
cumulative energy spectrum of a TGF sample of 130 
events, satisfying the same stringent selection crite- 
ria, detected during the period June 2008 - January 
2010. The relevant cumulative background as a func- 
tion of energy is calculated for events detected in the 
time interval T +l sec - T Q +21 sec, where T Q is the 
TGF start time. This method takes into account the 
20% orbital modulation of the MCAL background in 
a satisfactory way. We note that the MCAL back- 
ground level on TGF timescales is quite low (because 
of the MCAL anticoincidence vetoing, and of the satel- 
lite equatorial orbit), being 0.35 events/ms on aver- 
age. The MCAL detector response for different off- 
axis angles has been derived by a combination of sim- 
ulations and calibration data obtained up to a few 
MeV with radioactive sources, and up to 460 MeV 
at the beam test facility of the National Laborato- 
ries of Frascati (Italy). In-orbit calibration consis- 
tency checks were obtained for different Gamma-Ray 
Bursts at different angles. Remarkably, we find that 
the TGF spectrum extends up to 100 MeV with no 
exponential attenuation. Our data show the existence 
of a high-energy spectral component in addition to 
the well-known power-law (PL) component extending 
up to ~ 10 MeV. The additional component consti- 
tutes ~10% of the total emitted energy. A broken PL 
fit of the two components gives a differential photon 
energy flux F(E) - E -°- 5±0 1 for 1 MeV <E< E c , 
and F(E) ~ E^ 2J±0A for E c < E < 100 MeV, with 
E c = (7.1 ±0.5) MeV (all quoted parameter errors are 
1-a). 

Terrestrial Gamma- Ray Flashes turn out to be very 
efficient particle accelerators in our atmosphere. Our 
detected power-law emission between 10 MeV and 100 
MeV is difficult to reconcile with current RREA mod- 
els 0, 0, l20j - [23| . Some of these models are character- 
ized by acceleration over typical distances near, e.g., 
stepped-leader lightning sizes (~ 50 — 100 m) that 
correspond to a small number of avalanche lengths. 
On the contrary, an observed photon energy of 100 
MeV implies a lower limit on the acceleration distance 
d m i n ~ (1 km) (i?ioo) _1 where E\qq is the average 
electric field in units of 100 kV/m. These large-scale 
sizes are a significant fraction of the intra-cloud or 
cloud-to-ground distances over which potential drops 
of order of 100 MV can be established in thunder- 
storms. Furthermore, the detection of TGF emission 
in the 10-100 MeV range renews the interest for the 
neutron production in these energetic events as well 
as in normal lightning. A deteiled report on this dis- 
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FIG. 1: The background subtracted cumulative counts 
spectrum of the 130 TGFs detected by AGILE-MCAL 
during the period June, 2008 - January, 2010. The blue 
line shows the broken PL fit (see text), and the red 
curve is a pre-AGILE phenomenological model F(E) ~ 



E~ 



of index a = 0.4 ± 0.2, and exponential 



cutoff energy E c = 6.6 ± 1.2 MeV. Figure from 



covery is reported in [I'l 
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B. Gamma-ray localization of TGFs 

Since the discovery that TGFs can exhibit energies 
well above 20 MeV, we started searching for TGF 
events in the GRID, the AGILE gamma-ray imager 
sensitive above 20 MeV. The added value provided by 
the GRID is due to its imaging capabilities so that, 
for the first time, a TGF could be not only detected, 
but also localized from space directly in gamma-rays. 

In the period between June 2008 and Decem- 
ber 2009 the MCAL instrument triggered 119 bursts 
identified as TGFs according to the selection criteria 
discussed in [r|. For each of these bursts, the GRID 
dataset was searched for quasi-simultaneous gamma- 
ray events within a 200 ms time-window centered at 
the TGF start time To, defined as the time of the first 
MCAL-photon associated with the TGF. A peak in 
the cumulative distribution of the arrival times is ev- 
ident for the 2 ms time bin immediately following To. 
This peak includes 13 events, and the probability for 
it to be a statistical fluctuation (13 events or higher) 
is 6.5T0 -10 if we assume that GRID events are not 
correlated to TGFs and are distributed according to 
the Poisson law with the measured average rate of 
5.1 counts/s. All these GRID events take place dur- 
ing the TGF emission time interval estimated from 
MCAL data only. 

Figure [5] shows the scatterplot of the GRID events 
projection with respect to the AGILE footprint (the 
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FIG. 2: (a) scatterplot of the GRID events projection with 
respect to the AGILE footprint. The dashed circles are 300 
and 600 km in radius. The full AGILE visibility region has 
a 2600 km radius, (b) occurrence density vs. distance from 
footprint (each bin has been divided by the subtended area 
in km 2 ). Figure from [la ]. 



point at the Earth's surface on the straight line joining 
the satellite and the Earth's centre) , and the distribu- 
tion of the occurrence density vs. distance from foot- 
print (each bin has been divided by the subtended area 
in km ) . The GRID photon directions are obtained by 
the standard analysis pipeline for photons within the 
AGILE field of view, and by a custom version of the 
same software for photons outside the field of view. 
The photon directions are then back-projected from 
the satellite position, obtained by GPS data, to the 
Earth's surface, parameterized as the World Geode- 
tic System WGS-84 ellipsoid, in order to obtain the 
production site location. All 9 events are contained 
within a 1.14 sr solid angle, a factor 3.4 smaller than 
the solid angle subtended by the Earth at the satellite 
altitude of 540 km, which corresponds to a maximum 
visibility projected distance radius of ~ 2600 km from 
the satellite footprint. If the GRID events originate 
directly from the TGF production site, these results 
are consistent with a distance to footprint less than 
- 300 km estimated for RHESSI TGFs using sferics 
data [IB-Ill . 

The AGILE Tracker detections, described in details 
in [ll|, provide indeed the first direct TGF localiza- 



tion from space. The GRID detection of TGFs above 
20 MeV in about 7% of the cases confirms that TGFs 
emit a substantial amount of energy above this value, 
as reported in the previous section. 



C. TGF correlation with global lightning 
activity 

Correlation of TGFs with lightning activity was es- 
tablished shortly after their discovery. Several de- 
tailed one-to-one correlations of TGFs with individual 
lightning strokes have been published, see [HI, for 
some recent results, but a global climatology of TGFs 
is still missing, although some steps towards this di- 
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rection have been reported [lo|, HJ . 

Thanks to its detection capabilities and its very low 
inclination orbit, AGILE guarantees a high exposure 
above a narrow equatorial belt. Indeed, with about 
120 TGFs/year, the AGILE satellite record the high- 
est TGF detection rate surface density. This capabil- 
ity allowed a detailed high spatial resolution study of 
the correlation between TGFs and global lightning ac- 
tivity as detected from space by the Optical Transient 
Detector (OTD), onboard the Orbview-1 spacecraft, 
and the Lightning Imaging Sensor (LIS), onboard the 
TRMM NASA mission , (LIS/OTD public data avail- 
able at http://ghrc.msfc.nasa.gov/). The aim of this 
correlation study, based on mono and bi-dimensional 
Kolmogorov-Smirnov test, is to link TGFs to aver- 
age properties of lightning geographical distribution. 
Figure [SJ- A shows the annual lightning flash rate den- 
sity per km 2 per year with a spatial resolution of 
0.5 degrees for both longitude and latitude, based on 
LIS/OTD High Resolution Full Climatology (HRFC) 
data comprising about 10 years of observations (1995- 
2005) [13 ■ Figure [SB shows the MCAL exposure 
map, accounting also for the periods during which the 
trigger logic was active, while figure [3J-C shows the 
global lightning distribution multiplied by the MCAL 
exposure, that is directly comparable with AGILE 
TGF distribution. The bin size of 2.5 deg in longitude 
and 1 deg in latitude corresponds, on ground, to about 
275 km and 110 km respectively. In the same image all 
AGILE TGF positions are also shown, with white and 
black crosses. Finally, in Figure [5J-D the TGF (red) 
and the exposure-corrected lightning (black) longitude 
distributions are shown, summed over all latitudes. 

At a qualitative level, the longitudinal distributions 
for TGF and lightning above the continental areas 
(Figure [3J-D) show quite good agreement, and it is 
possible to recognize the main features of the three 
continental areas, like the sharp cut over the Congo or 
the double-peaked feature due to Sumatra and Borneo 
islands. However, to assess the compliance between 
the TGF and lightning distributions in a quantita- 
tive way, we used the Kolmogorov-Smirnov test (KS 
hereafter), which provides a probability value (P) for 
the null hypothesis that two unbinned data sets are 
drawn from the same distribution. After correction 
for the AGILE exposure, we find that we cannot con- 
sider the global TGF distribution as a random sample 
of the lightning distribution. Moreover, we find signif- 
icant regional differences in the degree of correlation. 
In particular, in the case of south east Asia we find a 
87% probability for the TGF distribution being a sub- 
sample of lightning, while this probability surprisingly 
drops to 3% in the case of Africa. 

The results shown here are an independent con- 
firmation of those reported in the work based on 
RHESSI data [13| concerning an excess of TGFs above 
central America and south east Asia and a correspond- 
ing depletion above Africa. However, reference [30| 



reports an eastward shift between TGF and light- 
ning distribution above south east Asia, which is cur- 
rently not supported by the the very good agreement 
we found for the two distributions above that region, 
but we must note that our result is obtained over a 
much narrower latitude region than that considered in 
that work. Moreover we must consider that the TGF 
dataset considered span very different periods of time. 
In fact, while AGILE dataset includes 118 TGFs de- 
tected in 12 months, the I s * RHESSI TGF catalog Q 
includes 144 TGFs in the same latitude range, but 
detected during 102 months of observation. So, while 
RHESSI TGFs are up to now the most complete sam- 
ple for climatological studies at a global scale, AG- 
ILE exhibits the higher detection rate surface density 
(TGFs/month per square degree) over a narrower ge- 
ographical region limited by its orbital inclination, al- 
lowing detailed studies of this climatic region. 

Based on the crude assumption that the observed 
TGF/lightning ratio holds at all latitudes we can also 
estimate a global rate of N ~ 220 570 TGFs per 
day, in agreement with previous estimates [H|. The 
correlation study presented in this section is described 
in details in reference [l9| . 

III. CONCLUSIONS 

AGILE is succesfully observing TGFs since June 
2008 and is currently one of the only three operating 
space instruments capable of detecting TGFs. The 
AGILE payload is very well suited for TGF science. 
Its main strength points can be summarized as follows: 

- AGILE-MCAL effective area peaks in the MeV 
range, the range where most of the TGF energy 
is radiated; 

- the MCAL energy range is extended up to 
100 MeV, allowing to probe the high energy tail 
of the TGF spectrum; 

- the trigger logic on time scales as short as 1ms 
and 290/is, well matching the TGF typical time 
scale, makes the AGILE sample not biased to- 
ward the brightest/longest events; 

- the MCAL design strategy, spatial segmen- 
tation in several independent detection units, 
makes the instrument less sensitive to dead-time 
and pile-up effects than monolithic detectors of 
equivalent volume; 

- event data with 2[is timing accuracy are avail- 
able for triggered events: time binning is limited 
by counting statistics only; 

- absolute timing accuracy better than 100/is al- 
lows precise timing for correlation with on- 
ground observation of sferic waves associated to 
lightnings; 
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FIG. 3: Lightning and TGF maps. A: LIS-OTD high resolution full climatology flashes rate [nashes/km 2 /year] (0.5x0.5 
deg per bin); B: MCAL exposure map [sec/bin] (2.5x1.0 deg per bin); C: LIS-ODT multiplied by MCAL exposure 
[flashes/km 2 ] (2.5x1.0 deg per bin). The crosses indicate the AGILE-TGF locations. The red borders indicate the 
continental zones considered in the analysis; D: Longitude distributions, summed over all latitudes, of the AGILE-TGF 
map (red) and LIS/OTD corrected map (black), corresponding to the maps showed in panel[3]-C, normalized to the total 
number of TGFs. Figure from [lgj ] 



- the AGILE-GRID trigger logic is sufficiently 
flexible to collect also high-energy photons com- 
ing from the Earth; 

- the AGILE orbit with 2.5° inclination is optimal 
for mapping the equatorial region, where most of 
the TGFs take place, with exposure much larger 
than other missions; 

Thanks to these capabilities, several important re- 
sults on TGF science based on AGILE observations 
have already been published [ilj - fl9l ], concerning the 
yet poorly understood TGF high-energy component, 
the first localization of TGFs from space, and the cor- 
relation of TGFs with global lightning activity. 
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